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Abstract. The derivation of the general solutions for stationary and static 
cylindrically symmetric Einstein spaces of Lewis form is revisited and the physical 
and geometrical meaning of the parameters appearing in the resulting solutions 
are investigated. It is shown that three of the parameters (and the value of the 
cosmological constant) are essential, of which one characterizes the local gravitational 
field and appears in the Cartan scalars, while the remaining two give information 
about the topological identification made to produce cylindrical symmetry. Other 
than the cosmological constant, they can be related to the parameters of the vacuum 
Weyl and Lewis classes, whose interpretation was previously investigated by da Silva 
C*** . ct al. (1995a, 1995b). 
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Introduction 



In a recent paper, Santos (1993) gave a set of solutions for stationary cylindrically 
symmetric spacetimes of the Lewis form (1932) with a cosmological constant. 
Krasinski had previously given these in a different form (Krasinski 1975a, 1994). Both 
use coordinates in which the static limit is hard to obtain. We shall reconsider the 
derivation of the complete set of stationary and static solutions. In another set of 
recent papers, da Silva et al. (1995a, 1995b) have given some physical identifications 
for the parameters of the corresponding vacuum solutions. Here we extend that 
investigation to the cases with cosmological constant. 

We first make some remarks on the problem of definition of cylindrical symmetry. 
In Kramer et al (1980) a spacetime is defined to be cylindrically symmetric if it is 
axisymmetric about an infinite axis and translationally symmetric along that axis. 
The definition of axisymmetry gives some difficulties (Mars and Scnovilla 1993, 1995). 
In the non-singular case, we can characterize axisymmetry by the requirement that 
some Killing vector (a) has closed spacelike orbits near the axis and (b) has length 
zero at the axis. In order for the axis to be non-singular we require elementary flatness 
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at the axis. If such an axisymmetric spacetime is also stationary, then one can prove 
that the stationary and axisymmetry Killing vectors commute and form surfaces with 
orthogonal transitivity, but these results depend on the existence of regular axis points 
(see chapter 17 of Kramer et al (1980)). 

These ideas are not currently well-defined if the axis is not regular, or has no 
regular part, but to exclude metrics without regular axes would be too restrictive as 
it excludes many of the known solutions which one would want to call cylindrical. In 
particular it rules out solutions which if continued to the axis (p = in the usual 
coordinate systems) would have strange behaviour but which can be joined at some 
P = Po > to a perfectly regular and well-behaved interior solution. This requirement, 
together with the existence of two commuting spacelike Killing vectors of which one 
has closed orbits and the other has non-closed orbits, so that they are transitive on 
surfaces of topology S 1 x R , could be taken as the definition of cylindrical symmetry. 
However, at present we know no way to test whether a metric could form a part of 
a globally well-behaved cylindrically symmetric solution in the sense of Kramer et al 
and it is therefore simpler to adopt the less demanding definition requiring only the 
existence of spacelike surfaces of symmetry with topology S 1 xR 1 . 

If there is no regular axis, one could define the presence of a singular axis by 
requiring the spacelike Killing vector with closed orbits to have a length which 
tends to zero, but the subsequent treatment of such spacetimes has some difficulties. 
In particular one cannot prove the orthogonal transitivity property, either in the 
stationary case (where it can occur with a timelike surface of transitivity), or 
in the more general non-stationary cylindrically symmetric case, in which Kramer 
et al (1980) assume it to occur with a spacelike surface of transitivity (see equation 
(20.1) there). In particular, if in the latter form one assumes that the metric 
coefficients are independent of time, one would wish to call the solution static 
even though the metric does not take the usual form for a stationary axisymmetric 
metric with orthogonal transitivity: for some such solutions see Chitre et al (1975), 
MacCallum (1983). Another possible problem is shown by the solution discussed by 
Lemos (1994) in which the "axis" p = in our coordinates is interpreted as the horizon 
of a cylindrical black hole. 

Thus the Lewis metric form used by Santos does not include all possible stationary 
Einstein metrics with cylindrical symmetry, defined as above. Even in the vacuum case 
the solutions of the Lewis form in general do not have well-defined axes (where the 
Killing vector with closed orbits has zero length), or if they do the axis is not in general 
elementarily flat. We therefore begin by simply assuming the Lewis form of the metric 
(Lewis, 1932), 

ds 2 = - fdt 2 + 2kdtdp> + tdip 2 + c u dr 2 + c^dz 2 (1.1) 

where /, k, £, v and /x are functions of r only, and the ranges of the coordinates t, r, 
z and p are initially taken to be 

— oo < i < oo — oo < r < oo — oo < z < oo < p < 2ir (1-2) 

with the hypersurfaces <p = and tp = 2n identified to ensure the cylindrical symmetry. 
Note that we have not assumed < r, in order both to avoid assuming the existence of 
an axis and to allow more freedom of coordinate choice. The coordinates are numbered 

x° = t x 1 =r x 2 = z x i = ip. (1.3) 

We should note immediately that the identification made on ip makes this an 
improper coordinate system: it does not obey the differential geometric requirement 
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of giving a homeomorphism between a region of the manifold and an open set in R . 
For that reason, we should really use at least two coordinate patches to cover the 
whole manifold. The permissible local transformations of coordinates, on these two 
or more patches, which are compatible with a identification in ip, with period P' say, 
are (MacCallum, 1997) 

t = Ti* + T , f = f(r), z = Z lZ + Z 0l <p = Xit + X 2 p + X , (1.4) 

where the upper case letters denote constants, and f is an arbitrary function. It is 
clear that the changes of origin of t, z and p do not alter any metric coefficient and 
can be ignored, so we need only consider the more restricted transformations 

t = T 1 t, r = f(r), z = Zz, <p = X 1 t + X 2 p. (1.5) 

However, we could also locally allow the more general form 

t — T\t + T 2 p, f = f(r), z — Zz, ip — X\t + X 2 ip. (1-6) 

In MacCallum (1997), it is shown that of the four parameters I\, T 2 , X\ and X 2 
only two are essential; they fix the topological identification needed to reach the 
actual solution starting from its locally equivalent standard form. If in the above 
transformations the barred coordinates refer to this standard form then the essential 
parameters S and P of MacCallum (1997) are —X 2 /T 2 and 2nX 2 (or we could say, 
more succinctly, that X 2 and T 2 are essential). 

The radial coordinate can be chosen so that v = p, (as Santos (1993) does initially) 
or so that v = 0, which is the more convenient choice of radial coordinate introduced 
later by Santos (1993), called r* in that paper, which is in fact the natural r coordinate 
arising from the existence of hypersurface-homogencity of Bianchi type I on timelike 
surfaces (which thus necessarily have geodesic normal unit vectors). With this choice 
the remaining coordinate freedom of r is f — ±r + R where R is a constant. 

2. The Einstein space solutions 

We now revisit and complete the derivation of the solutions of the Einstein equations 

Rab = A. 9ah (2.1) 

for the metric form considered. Following Santos (1993) we find that if u = f/£ and 
v = k/£ then with the choice of r given by v = fi, and denoting d/dr by ', we have 

u' = &{u + v 2 ), v' = <£'(u + v 2 ) (2.2) 

where (pQ 1 )' = = (p&'Y , p (denoted D in Santos (1993)) being the usual radius 
parameter p 2 = f£ + k 2 (which in vacuum and some other cases can also be taken as 
an isothermal coordinate in the (r, z) surfaces). There are now several possibilities. 
If both 0' and $' are non-zero, they must be multiples of one another, say $' = AO', 
and this leads (as in Santos (1993)) to v = Au + B where A and B are constants, i.e. 
to k = Af + B£. If one of 0' and <£' is zero, so that either O or $ is constant (possibly 
zero), then again /, k and £ are linearly dependent, and this is also true if £ = so that 
u and v are not well-defined. Thus in all cases the proposition of MacCallum (1997) 
applies, and all solutions are locally equivalent either to a diagonal metric, or one 
with a null Killing vector, or a 'windmill ' form ( which can be considered as a complex 
diagonal metric). Moreover, the solution for the diagonal case will give the solutions 
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for the remaining cases if we allow complex values for parameters and coordinates and 
take limits as in Mcintosh (1992). 

Taking the form (1.1) with k = = v wc find (cf. Linet (1985) and Santos (1993)) 
that if G = pe^ 2 then 

G" + SAG = 0. 

The distinct solutions (after choice of the origin and sign of r) are 



G = C cosh (^^/3\A\r^ , A < 0; (2.3) 

= Gexp (v / 3|Afr) , A < 0; (2.4) 

= Gsinh (vT|A[r) , A < 0; (2.5) 

G = Cr, A = 0; (2.6) 

G = C sin (V3Ar) , A > 0, (2.7) 



where in all cases C > 0. The interpretation of G < is unclear. For example, it 
might represent a double-covering of the spacetime or a region interior to a horizon. 
The vacuum solutions are included for comparison purposes although these solutions, 
which can all be related to the Kasner form (see Mcintosh (1992) and references 
therein), have been discussed in many earlier papers. By scaling the three ignorable 
coordinates we can arrive at C = l/^/|3A| if A ^ 0: for A = wc can normalize so 
that C = 1, and we denote these normalized cases by Go- They have the property 
that 

(G ) 2 + 3A(G ) 2 = r, 

where in the five cases (2.3-2.7) r\ is respectively — 1, 0, 1, 1 and 1. 

The case A < gives three forms of solution, but Linet (1985) considers only (2.5), 
and the differences are not obvious in Krasinski's complicated form. In Santos (1993) 
all three are included in the general form 

G = Gi cosh (\/^3A>) + G 2 sinh (\/^3Ar) , (2.8) 

the three cases corresponding to \C\\ > IG2I, \C\ \ — | C2 1 , and \C\\ < |G 2 |. In the first 
and third cases the G above is \J\C\ — Gf|, while in the second case G = \C\\. The 
origin shifts in r required to reach (2.3) and (2.5) are respectively arctanh(G2/Gi) and 
arctanh(Gi/G2). In the similar formulae for the case (2.6) G = G2 and the shift is 
G1/G2 and in the case (2.7) G = y ' G\ + Gf and R = arctan(Gi/G 2 ). 

With these specializations, and using allowable scalings of the remaining 
coordinates, we can write the resulting metric (Linet, 1985) as 

ds 2 = dr 2 + Go /3 [ £ i ex P( 2 fe' - 1 /3)C/)(dx J ) 2 ] (2.9) 

J=0,2,3 

where Sj is 1 except for x° = t, for which it is —1, and dll/dr — 1/Go- (No constant 
of integration is given in U as it can be absorbed by rescalings of coordinates, but the 
rescalings allowed appear explicitly in the discussion below.) The constants pj (which 
are 1/3 + KKj/2 in terms of Linet's notation) obey 

I> = 1 ' E^ = (^ +1 )/ 3 - ( 2 - 10 ) 

3 3 
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The parameterization here has been chosen to agree with the usual forms of the 
Kasner vacuum metric (equations (11.50) and (11.51) in Kramer et al (1980)) and its 
'windmill' counterpart (Mcintosh, 1992) for the case A = 0, where G = r = cxp U. As 
in the Kasner vacuum case, a complex coordinate change gives related 'cosmological' 
solutions, i.e. solutions for which the essential variable has a timelike character. These 
related metrics have been given several times: Kasner himself (1925) gave the case 
related to (2.7), Saunders (1967) (see equation (11.52) in Kramer et al (1980)) the 
cases related to (2.5) and (2.7), Kale (1970) the case related to (2.3), Spindel (1979) 
gave the full set, and there may be other occurrences we have not found. The metric 
forms are included in those considered by Kellncr (1975), but he seems not to give 
explicitly the solutions for a cosmological constant. Linet (1985) only considered the 
forms with (2.5) and (2.7), as he wished to have an axis p — at some r. 

Note that for rj = the only real solution is pi = 1/3 for all i: other cases with 
rj = 0, and all cases with n = — 1, will have complex values of p^. It is not obvious 
whether in general the resulting metrics, with the coordinates also considered to be 
complex, will have a real section, as happens in the 'windmill' case (Mcintosh, 1992). 
However, there are clearly cases where this does happen: for example, with r/ = — 1, 
if po = 1 , pi = p~3 = V2i I y/3 where we get a real solution by replacing z and <p by 
x — iy and x + iy respectively. Linet ignores these cases, taking his K 2 > 0. They 
are implicitly present in Santos (1993), where they correspond to complex values of 5 
and cry//? in his equation (51), as corrected (see the erratum to the paper and below). 
The similar examples with p2 = 1 are covered by his equation (57). If we assume that 
P2 is real and po = a + ib = p~s, we find this is possible, with non-zero real b, for values 
of a outside [0, 2/3] if rj = 1, for a ^ 1/3 if n = and for any a if 7/ = —1. The 
cases with a null Killing vector can be obtained by going to coordinates in which the 
metrics appear stationary and then taking a suitable limit. 

There are clearly two essential parameters in (2.9), one of the p,, say P2, and 
the cosmological constant (and the discrete parameter rj if A < 0). One could argue 
that A should not be regarded as a parameter, since in the usual Einstein theory it 
would be a universal constant. However, if the term is regarded as arising as a relic of 
quantum processes, or if we consider that in writing down the whole family we have 
an infinite number of completely disjoint spacetimes, and hence no physical reason for 
using the same A, it may make sense to treat A as a parameter. To make the solution 
cylindrical, ip must be periodic, and the period is specified by another essential global 
parameter (see MacCallum (1997)): note that scaling ip to make the period 27r alters 
the g vv of (2.9) by a constant factor. Thus in the form (2.9) for a cylindrical metric 
with a periodic ip there are three parameters, A, p 2 and the period (plus the discrete 
parameter 77 if A < 0). 

For the distinct possible Go given by (2.3-2.7) the formulae for U are respectively 




A < 



A < 



A < 



(2.13) 
(2.14) 
(2.15) 



(2.11) 



(2.12) 



U 



lnr, 



A = 




A > 



cf. Santos (1993), equations (45)-(49). 
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We can now obtain the general forms of the corresponding classes of stationary 
cylindrically symmetric metrics by applying the transformations (1.6), taking the 
barred coordinates to be the ones for the standard form just discussed. This will 
introduce a fourth essential parameter (not counting rj) which specifies the lines in the 
original ((p, t) surface parallel to t = along which the topological identification is to 
be made, as well as some inessential parameters: as noted above X 2 and T 2 can be 
considered the essential parameters in addition to A and p 2 . We do the calculations 
explicitly only for the diagonal form, assuming p 2 to be real, which leaves us the 
possibility of complex po and p 3 . Then 

ds 2 = dr 2 + G 2 / 3 (Z 2 cxp(2(p 2 - l/3)U)dz 2 

+ \X 2 2 exp(2(p 3 - 1/3)17) - T| exp(2(p - l/3)*7)]d^ 2 (2.16) 
+ 2[X,X 2 cxp(2(p 3 - 1/3)17) - T{F 2 exp(2(p - l/3)C/)]didy> 
-[T 2 exp(2(p - 1/3)17) - X 2 cxp(2(p 3 - l/3)[/)]dt 2 ) 

and the full coordinate freedom includes replacing r in the formulae for Go and U by 
±r + R where R is a constant. The geometry has 4 essential parameters A, p 2 , X 2 
and T 2 . G is now given by 

G 2 = Z 2 {T X X 2 - T 2 X X ) 2 G 2 . 

The remaining free parameters T\, X\, R and Z are inessential in specifying the 
geometry of the region covered, but may be required in order to match this metric 
to an interior using the Lichnerowicz conditions, which are the same as the earlier 
Darmois conditions but expressed in 'admissible coordinates' continuous across the 
boundary (see Bonnor and Vickers (1981) and applications in Krasinski (1975a) and 
Bonnor et al (1997)). This coordinate choice, and therefore the inessential coordinate 
parameters, are not required in the Darmois form, in which the conditions are just that 
the first and second fundamental forms of the surface, considered from the two sides, 
must agree, and thus only invariantly-defined parameters can appear. Of course, on 
trying to match a given solution to the solution here one may find it has more essential 
parameters than can in fact be matched by the solutions here, or parameters which 
appear in an incompatible way: there is no guarantee that a match is possible. 

These solutions can now be identified with earlier forms. Krasinski (1975a) used a 
"radial" coordinate x 2 = k/f in the above notation: this leads to rather complicated 
forms, and the static limit k = is awkward to handle. (This coordinate choice is 
made in order to fit the exteriors easily to his non-static interior solutions.) Indeed 
Krasinski (1975b) gives only the case A = in the explicitly static form, as his equation 
(7.2), and says "The generalization of (7.2) to the case A 7^ is unexpectedly very 
involved, so I do not present it here" , though in fact the solutions themselves appear in 
stationary form. We can identify his classes, simply by looking at the coefficients in the 
linear dependence of /, k and I (cf. MacCallum (1997)) as follows. His class A is locally 
equivalent to the 'windmill' form, Class B is locally equivalent to a real diagonal form, 
and Classes C and D are locally equivalent to the form with a null Killing vector, the 
difference being that in Class D the identification is made along this Killing vector's 
orbit. The identification of Class B with the (locally) static cases other than his (9.3) 
is stated in Krasinski (1975a), in the paragraph following Table IX, and the proof of 
this (which was known to Krasinski at that time (private communication, 1997)) is 
hinted at but not completely described in the paragraph following his equation (9.2). 
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The version given originally by Santos (1993) contains some typographical errors, 
corrected in the subsequent erratum, 
(i) If K = -a 2 < 0, we have 

- = -sinh — , (2.17) 



P a \ 2/3 J 

k , faO\ (3 . , fa@\ 

= - cosh — - - sinh — , (2.18) 



p V 2/9 / a V 2/5 

/ „„ , /a9\ (3 2 + a 2 . , /a0\ 



(h) ifjr = o, 

£ _ _e 



(2.20) 



i~i-f, 

- = -£0-2/3. (2.22) 
P 2 



(hi) If X = a 2 > 0, 



^ 1 • ,„„„, 
- = - sin — (2.23) 

cos -r sin > ( 2 - 24 ) 



2/? cos — - sin — , (2.25) 



where 



p V 2(3 J a \ 2(3 
f faQ\ p 2 -a 2 
-P=- 2PC0S {W) —~\20J 

6 = 7 + C, e 3 ^ 2 =eGcxp((50), il = jdr/G. 



5, 7, e and ( are constants of integration. 

The 9 constant parameters, A, C\, C2, a, (3, 7, S, e and (, in these solutions are 
related, in each case, by one equation, as follows: 

(i) K < 



3|A|(C 2 -C 2 ) + -5 2 + _^ = 0, A<0; (2.26) 
-Cl + h 2 + ^-=0, A = Q . (2 27) 



-3A(C 2 + C 2 ) + I 5 2 + ^ = 0, A>0. (2.28) 

(ii) K = 

3 |A| (C 2 - C|) + J<5 2 - 0, A < 0; (2.29) 

C* 2 ±i<5= 0, A = 0; (2.30) 

-3A (C? + C 2 2 ) + J(5 2 = 0, A > 0. (2.31) 
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(iii) K > 

3 |A| (C x 2 - C|) +I^_3g^ =0] A<0; (2 32) 

-3A(Cf + C?) + ^ 2 -^gf = 0, A>0. (2.34) 

We now seek to identify the 9 parameters in the form of Santos (1993) with the 
parameters of the general solution (2.16). We have already identified C\ and C 2 with 
C and R above and we see that ft = f//3|A|C if A ^ 0, or f//C if A = 0. For the five 
possible G we have Z 2 (T 1 X 2 - T 2 Xi) 2 equal to, respectively 3|A|(C? - C|), 3|A|Cf , 
3|A|(C 2 - C 2 ), C 2 , and 3A(C 2 + C|). We find <5 = v /3|A|C , (3p 2 - 1) if A ^ or 
5 = C(3p 2 -l) if A = 0, and e = Z 2 /(TiX 2 - T 2 Xi). Incase (i) a 7 = 6|A|C/3(p 3 ~Po) 
(or erf = 2C(i(pz — Po) if A = 0); in case (iii) p an d Pz are complex conjugates and 
«7 = -6i|A|C/3(p 3 -p )- The relations listed as (2.26-2.28) and (2.32-2.34) are then 
equivalent to the second of the equations (2.10), the first of these being already built 
into the forms above. 

The parameters a and (3 are most easily identified from the linear dependence of 
/, k and I which for case (i) is 

/ - 2(3k - f3 2 £ = -a 2 l 

and can be compared with 

/ + # + ^-^ = ° 

giving us that j3 + a and (3 — a are T\ /T 2 and X\ / X 2 (the correct pairing depending 
on the choice which keeps t timelike). Finally, use of the transformation 

t = t- T 2 ^/T! , (p = -X 1 t/X 2 + if, (2.35) 

on the metric (2.16), and the same in terms of a and (3 on that of Santos, enables us 
to identify ( in terms of \n{X 2 /Ti). 



3. Invariant properties 

As discussed above, we have 4 essential parameters in the general stationary metric, 
including A but excluding r\. From the work of MacCallum (1997) and the results of 
Da Silva et al (1995a), Da Silva et al (1995b) we can expect that 2 of these (T 2 and 
X 2 ) will appear in the holonomy around the curves on which only tp changes, and the 
remaining two (A and p 2 ) in invariants of the local curvature. Further parameters 
may appear from coordinate conditions on a matching to an interior solution. 

We can compute the "Cartan scalars" which completely characterise the local 
geometry according to the procedure for invariant classification described in Karlhcde 
(1980), Paiva et al (1993) and MacCallum and Skea (1994). For case (i) with the G 
of equation (2.5) and C = l/y / 3|A|, this calculation, which was carried out using the 
computer algebra package CLASSI built on SHEEP, gave, as the non-trivial classifying 
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invariants, the cosmological constant and 

*o = *4 

= l/4a/rV|A|sinh- 2 (V3|A>) 

+ l/4a/3-S|A|cosh(v / ^r)sinh- 2 ( v /3jAfr-) (3.1) 
* 2 = -l/8a 2 /3- 2 7 2 |A| smh- 2 (^3\A\r) + 1/3|A| sinh- 2 ( y^r) 
- 1/65|A| cosh(y3|A|r) sinh- 2 (V3|A|r) 

O* l' = £>*50' 

= V6 (l/16a 3 /3- 3 7 3 |A| 3 / 2 sinh- 3 (v/3lA>) 

- 5/12a/T V|A| 3/2 cosh( y/3\A\r) sinlT 3 (v / 3iA|r) (3.2) 

- 5/24a/3~ 1 7 |A| 3 / 2 sinh-^VajAjr) 
-2/3a/3- 1 7 |A| 3/2 sinh- 3 (v^]A7r)) 

£>*io' = .0*41' 

= V6(-l/16a 3 /?- 3 7 3 |A| 3 / 2 sinh- 3 (v^v|r) 

- l/12a/T V|A| 3/2 cosh(y3|A|r) sinh- 3 (y3|A|r) (3.3) 

- l/24a/?- 1 7 |A| 3 / 2 Binh-^V^jAfr) 
+l/6a/r 1 7 |A| 3/2 sinh- 3 ( v ^]A7r)) 

0*21' = ^0*30' 

= V6 (l/8a 2 /?- 2 7 2 |A| 3 / 2 cosh(73lAJr) sinh- 3 ( v /3Rr) 

+ l/12(5|A| 3 / 2 sinh- 1 ( v /3]A]r) (3.4) 
+ l/65|A| 3 / 2 sinh- 3 ( v /3|A|r) 
-l/3|A| 3 / 2 cosh(y3lAJr) sinh" 3 (y/3\K\r)) 

We see that the other essential parameter that appears is cry//3 which is a linear 
combination of the pi and so equivalent to specifiying pi as expected. The results for 
the other possible cases described above are similar in form. Thus, as expected, we 
find that only two parameters are essential locally: all the other parameters can be 
locally removed by coordinate transformations. 

We can also relate the Cartan scalars above to their limiting forms for the Weyl 
vacuum class considered by da Silva et al. (1995b). In the limit A = we can 
identify the other parameters in the Cartan scalars with the n of the Weyl form 
by 5 = 2(n 2 - 3)/(n 2 + 3), cry/ (3 = 8n/(n 2 + 3), and with these relations the Cartan 
scalars above have the da Silva et al. values as their limits. 
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